SYMMETRIC o-STABLE SUBORDINATORS AND CAUCHY 

PROBLEMS 



ERKAN NANE 



Abstract. We survey the results in Nane (E. Nane, Higher order PDE's and 
iterated processes, Trans. American Math. Soc. (to appear)) and Bacumer, Meer- 
schaert, and Nane (B. Baeumer, M.M. Meerschaert and E. Nane, Brownian subor- 
dinators and fractional Cauchy problems: Submitted (2007)) which deal with PDE 
connection of some iterated processes, and obtain a new probabilistic proof of the 
equivalence of the higher order PDE's and fractional in time PDE's. 



1. Introduction 

In recent years, starting with the articles of Burdzy [9~| [TO], researchers have shown 
interest in iterated processes in which one changes the time parameter with one- 
dimensional Brownian motion. 

To define iterated Brownian motion Z t , due to Burdzy [9], started at z G K, let 
X t + , Xf and Y t be three independent one-dimensional Brownian motions, all started 
at 0. Two-sided Brownian motion is defined to be 



xt, t > 

Xf_ t)i t<0. 



Then iterated Brownian motion started at z G M is Z t = z + X(Y t ), t > 0. 

1.1. BM versus IBM. This process is not Markovian or Gaussian but, it has many 
properties analogous to those of Brownian motion. We list a few; 

IBM Z t has stationary (but not independent) increments, and is a self-similar 
process of index 1/4. Laws of the iterated logarithm (LIL) holds: usual LIL 
by Burdzy [9] shows 

Z(t) 2 5 / 4 

limSUp : ,, ; . = ^777 d.S. 

Chung-type LIL is obtained by Khoshnevisan and Lewis [19] and Hu et al. [TT] . 
Khoshnevisan and Lewis [18] extended results of Burdzy [10], to develop a stochastic 
calculus for iterated Brownian motion. Burdzy and Khosnevisan [12] showed that 
IBM can be used to model diffusion in a crack. Local times of this process was studied 



1991 Mathematics Subject Classification. 60J65, 60K99. 

Key words and phrases. Iterated Brownian motion, PDE connection, a— stable process, a-time 
process, Levy process, Brownian subordinator, Caputo derivative, fractional derivative in time. 



by Burdzy and Khosnevisan [TT], Csaki, Csorgo, Foldes, and Revesz [13], Shi and Yor 
|32j . Xiao [51], and Hu [TB]. Baiiuelos and DeBlassie [B] studied the distribution 
of exit place for iterated Brownian motion in cones. DeBlassie [H] studied the 
lifetime asymptotics of iterated Brownian motion in cones and Bounded domains. 
Nane EH EHl EHl G9] extended some of the results of DeBlassie. 



1.2. PDE connection. In addition to the above properties of IBM, there is an inter- 
esting connection between iterated Brownian motion and the biharmonic operator 
A 2 . Allouba and Zheng [2] show that if we replace the outer process X(t) in the 
definition of iterated Brownian motion with a continuous Markov process with L x 
as the semigroup generator, then u(t,x) = E x [f(Z t )] := E[f(Z t )\Z = x] solves the 
Cauchy initial value problem 

(1.1) ^-u(t,x) = La 'ffl + L x 2 u(t,x); u(0,x) = f(x), t > 0,x G K d . 



When Z t is an iterated Brownian motion, this was also obtained by DeBlassie [H] by 
a different method. Let Z\ = X(\Y t \), then Allouba and Zheng [I] shows E x [f(Z t )] = 
E x [f(Zl)). 



1.3. Fractional Cauchy problems. Nigmatullin [30] gave a physical derivation of 
the fractional kinetic equation for some special (3 

(1.2) -^u(t,x) = L x u(t,x); u(0,x) = f(x) 

where < (3 < 1 and L x is the generator of some continuous Markov process X (t) 
started at x = 0. Here d^g(t)/dt" is the Caputo fractional derivative in time, which 
can be defined as the inverse Laplace transform of s^g^s) — s /3_1 (?(0), with g(s) = 
Jo°° e ~ st g(t)dt the usual Laplace transform. 

Mathematical study of equation (II. 2p was initiated by [321 EDI EI] • The existence 
and uniqueness of solutions to equation (jl.2p was proved in [201 EI] • This equation 
was also used by Zaslavsky [35] for Hamiltonian chaos. 

Baeumer and Meerschaert [1] and Meerschaert and Scheffler [23J show that the frac- 
tional Cauchy problem (11.21) is related to a certain class of subordinated stochastic 
processes; take D t to be the stable subordinator, a Levy process with strictly increas- 
ing sample paths such that E[e~ sDt ] = e~ tsP , see for example Bertoin [Sj. Define the 
inverse or hitting time or first passage time process 

(1.3) E t = inf{x > : D(x) > t}. 

The subordinated process Z t = Xo(E t ) occurs as the scaling limit of a continuous time 
random walk (also called a renewal reward process), in which iid random jumps are 
separated by iid positive waiting times (Meerschaert and Scheffler (2004) [23]). The- 
orem 3.1 in Baeumer and Meerschaert [1] shows that, in the case p(t,x) = T(t)f(x) 



is a bounded continuous semigroup on a Banach space, the formula 

/OO £ POO £ 

p((t/sf, x)gp(s) ds = - J p(x, s)g l3 {^ T jp)s- 1/(3 - 1 ds 

yields a solution to the fractional Cauchy problem fll.2p . Here, gp(t) is the smooth 
density of the stable subordinator, with gp{s) = J °° e~ st gp{t) dt = e~ sl3 . 

2. Brownian subordiantors and fractional Cauchy problems 



We give a probabilistic proof of the following theorem. A variation of this result 
was realized by Orsingher and Benghin [21] for a version of iterated Brownian motion. 

Theorem 2.1 (Baeumer, Meerschaert and Nane (2007) [5]). Let L x be the generator 
of a Markov semigroup T(t)f(x) = E x [f(X t )}, and take f G D(L X ) the domain of 
the generator. Then, both the higher order Cauchy problem U.l\) and the fractional 
Cauchy problem U.Sfy with (3 = 1/2, have the same solution 

(2.1) u(t,x) = E x {f{Z t )] = -j=J^ T(s)f(x)e W (--)ds. 

Proof. E t is the inverse of a 1 — 1/a stable subordinator. E t then is the local time 
of symmetric stable process of index a. In the case a = 2, local time of Brownian 
motion is the same as sup 0<s<t B s . On the other hand, sup 0<s<t 5 s and \B t \ are 
same in distribution by the reflection principle. Hence, E t and \B t \ have the same 
one-dimensional distributions, implying the result of the theorem. □ 

We obtain the following corollary of our theorem 

Corollary 2.2 (Baeumer, Meerschaert and Nane (2007) [5]). For any continuous 
Markov process X(t), both the Brownian-time subordinated process X(\Y t \) and the 
process X(E t ) subordinated to the inverse 1/2-stable subordinator have the same one- 
dimensional distributions. Hence they are both stochastic solutions to the fractional 
Cauchy problem U.fy) . or equivalently, to the higher order Cauchy problem (II. ip . 

In contrast to the previous corollary, we have 

Theorem 2.3. Let Y be a symmetric stable process of index 1 < a < 2, and E t is the 
inverse of a stable subordinator of index 1 — 1/a. The processes X(E t ) and X(\Y t \) 
do not have same one- dimensional distribution. 

Proof. Let V\ be the local time at x = 0. L\ has the same one-dimenional distributions 
as E t . Lemma 1 in Hawkes [15] implies that 

(2.2) P[L\ > X] ~ d\- a/2 exp(-C ah X a ). 
Proposition 4 in Bertoin [S] shows 

(2.3) P[Yi > u] ~ P[ sup Y s > u] ~ cu~ a . 

0<s<l 



The results in equations (12. 2p and (12.31) establish that in the case Y t is a symmetric 
stable process of index a < 2, \Y t \ and E t do not have the same one-dimensional 
distributions. □ 

When the outer proces is Levy process we have uniqueness of the solutions in 
Theorem 12. 1[ The proof relies on a Laplace- Fourier transform argument. 

Theorem 2.4 (Baeumer, Meerschaert and Nane (2007) [5]). Suppose that X(t) = 
x + Xo(t) where Xo(t) is a Levy process starting at zero. If L x is the generator of 
the semigroup T(t)f(x) = E x [(f(X t ))] on L 1 {R d ), then for any f G D(L X ), both the 
initial value problem (II. ip . and the fractional Cauchy problem (ll.2p with (3 = 1/2, 
have the same unique solution given by (12.11) . 

An easy extension of the argument for Theorem 12.41 shows that, under the same 
conditions, for any n = 2, 3, 4, . . . both the Cauchy problem 

(2-4) — ^ = £ f(jM Li/(x) + L>( *' X)] M(0 ' X) = f(x) 

and the fractional Cauchy problem ( j 1.2ft with (3 = 1/n have the same unique solution 
given by u(t,x) = J °° p((t/s) l3 ,x)gp(s) ds with (3 = 1/n. Hence the process Z t = 
X(E t ) is also the stochastic solution to this higher order Cauchy problem. 

3. Other subordinators 

An a-time process is a Markov process subordinated to the absolute value of an 
independent one-dimensional symmetric a-stable process: Z t = B(\S t \), where B t is 
a Markov process and S t is an independent symmetric a-stable process both started 
at 0. Let Z% = x + Z t the process started at x. 

This process is self similar with index l/2a when the outer process X is a Brownian 
motion. In this case, Nane [27] defined the local time of this process and obtained 
laws of the iterated logarithm for the local time for large time. 

3.1. PDE-connection: 

Theorem 3.1 (Nane (2005) [26]). Let T(s)f(x) = E[f(X x (s))] be the semigroup of 
the continuous Markov process X x (t) and let L x be its generator. Let a = 1. Let f be 
a bounded measurable function in the domain of L x , with Dijf bounded and Holder 
continuous for all 1 < i, j < n. Then u(t,x) = E[f(Zf)\ solves 

^-ru(t,x) = - 2Lx fj x ^ _ L 2 x u(t,x); u(0,x)= f(x). 
at 1 nt 

For a = l/m^l rational: the PDE is more complicated since kernels of symmetric 
a-stable processes satisfy a higher order PDE: 



We also have to assume that we can take the operator out of the integral. This is 
valid for a — 1/m, m — 2, 3, • • * , by a Lemma in Nane |26j . 

Theorem 3.2 (Nane (2005)[2S])- Let a 6 (0,2) be a rational with a = l/m, where 
I and m are relatively prime. Let T(s)f(x) = E[f(X x (s))} be the semigroup of the 
continuous Markov process X x (t) and let L x be its generator. Let f be a bounded 
measurable function in the domain of L x; with D 1 f bounded and Holder continuous 
for all multi index 7 such that \^\ = 21. Then u(t,x) = E[f(Zf)] solves 

f) 2m 1 / f)2l-2i \ 

(-!) m ^M = " 2 E ^ (0,S)U ) L T l m - Ll l u(t,x); 
w(0,a;) = f(x). 

For some other connections of PDE's and iterated processes, see papers by Nane 
and Allouba and Zheng [2], Allouba [I], Baeumer et al. [S] and references therein. 

4. Open Problems 

Question 1. Looking at the governing PDE for subordinators other than Brownian 
motion, are there any fractional in time PDE which has the same solution as the 
higher order PDE? 

Question 2. Are there PDE connections of the iterated processes in bounded domain 
as the PDE connection of Brownian motion in bounded domains? 
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